Anthrax toxin lethal factor domain 3 is highly mobile and responsive to ligand binding.

The secreted anthrax toxin consists of three components: the protective antigen (PA), edema factor (EF) and lethal factor (LF). LF, a zinc metalloproteinase, compromises the host immune system primarily by targeting mitogen-activated protein kinase kinases in macrophages. Peptide substrates and small-molecule inhibitors bind LF in the space between domains 3 and 4 of the hydrolase. Domain 3 is attached on a hinge to domain 2 via residues Ile300 and Pro385, and can move through an angular arc of greater than 35° in response to the binding of different ligands. Here, multiple LF structures including five new complexes with co-crystallized inhibitors are compared and three frequently populated LF conformational states termed `bioactive', `open' and `tight' are identified. The bioactive position is observed with large substrate peptides and leaves all peptide-recognition subsites open and accessible. The tight state is seen in unliganded and small-molecule complex structures. In this state, domain 3 is clamped over certain substrate subsites, blocking access. The open position appears to be an intermediate state between these extremes and is observed owing to steric constraints imposed by specific bound ligands. The tight conformation may be the lowest-energy conformation among the reported structures, as it is the position observed with no bound ligand, while the open and bioactive conformations are likely to be ligand-induced.

Analysis of the Errors in the Electrostatically Embedded Many-Body Expansion of the Energy and the Correlation Energy for Zn and Cd Coordination Complexes with Five and Six Ligands and Use of the Analysis to Develop a Generally Successful Fragmentation Strategy.

In the present paper, we apply the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE) to the calculation of zinc-ligand and cadmium-ligand bond dissociation energies, and we analyze the errors due to various fragmentation schemes in a variety of neutral, positively charged, and negatively charged Zn2+ and Cd2+ coordination complexes. As a result of the analysis, we are able to present a new, simple, and unambiguous fragmentation strategy. Following this strategy, we show that both methods perform well for zinc-ligand and cadmium-ligand bond dissociation energies for all systems studied in the paper, including a model of the catalytic site of the zinc-bearing anthrax toxin lethal factor (LF), which has garnered substantial attention as a target for drug development. To draw general conclusions we consider ten pentacoordinate and hexacoordinate zinc and cadmium containing coordination complexes, each with 10 or 15 different fragmentation schemes. By analyzing errors, we developed a prescription for the optimal fragmentation strategy. With this scheme, and using MP2 correlation energies as a test, we find that the electrostatically embedded three-body expansion of the correlation energy (EE-3B-CE) method is able to reproduce all 53 conventionally calculated bond energies with an average absolute error of only 0.59 kcal/mol. The paper also presents EE-MB-CE calculations using the CCSD(T) level of theory on an LF model system. With CCSD(T), EE-3B-CE has an average error of 0.30 kcal/mol.

Identification of novel anthrax toxin countermeasures using in silico methods.

Anthrax is an acute infectious disease caused by the spore-forming, gram-positive, rod-shaped bacterium Bacillus anthracis. The anthrax toxin lethal factor (LF) is the primary anthrax toxin component responsible for cytotoxicity and host death and has been a heavily researched target for design of postexposure therapeutics in the event of a bioterror attack. Various computer-aided drug design methodologies have proven useful for pinpointing new antianthrax drug scaffolds, optimizing existing leads and probes, and elucidating key mechanisms of action. We present a selection of in silico virtual screening protocols incorporating docking and scoring, shape-based searching, and pharmacophore mapping techniques to identify and prioritize small molecules with potential biological activity against LF. We also recommend screening parameters that have been shown to increase the accuracy and reliability of these computational results.

Development of a comprehensive, validated pharmacophore hypothesis for anthrax toxin lethal factor (LF) inhibitors using genetic algorithms, Pareto scoring, and structural biology.

Anthrax is an acute infectious disease caused by the spore-forming bacterium Bacillus anthracis. The anthrax toxin lethal factor (LF), an 89-kDa zinc hydrolase secreted by the bacilli, is the toxin component chiefly responsible for pathogenesis and has been a popular target for rational and structure-based drug design. Although hundreds of small-molecule compounds have been designed to target the LF active site, relatively few reported inhibitors have exhibited activity in cell-based assays, and no LF inhibitor is currently available to treat or prevent anthrax. This study presents a new pharmacophore map assembly, validated by experiment, designed to rapidly identify and prioritize promising LF inhibitor scaffolds from virtual compound libraries. The new hypothesis incorporates structural information from all five available LF enzyme-inhibitor complexes deposited in the Protein Data Bank (PDB) and is the first LF pharmacophore map reported to date that includes features representing interactions involving all three key subsites of the LF catalytic binding region. In a wide-ranging validation study on all 546 compounds for which published LF biological activity data exist, this model displayed strong selectivity toward nanomolar-level LF inhibitors, successfully identifying 72.1% of existing nanomolar-level compounds in an unbiased test set, while rejecting 100% of weakly active (>100 µM) compounds. In addition to its capabilities as a database searching tool, this comprehensive model points to a number of key design principles and previously unidentified ligand-receptor interactions that are likely to influence compound potency.

Electrostatically Embedded Many-Body Expansion for Neutral and Charged Metalloenzyme Model Systems.

The electrostatically embedded many-body (EE-MB) method has proven accurate for calculating cohesive and conformational energies in clusters, and it has recently been extended to obtain bond dissociation energies for metal-ligand bonds in positively charged inorganic coordination complexes. In the present paper, we present four key guidelines that maximize the accuracy and efficiency of EE-MB calculations for metal centers. Then, following these guidelines, we show that the EE-MB method can also perform well for bond dissociation energies in a variety of neutral and negatively charged inorganic coordination systems representing metalloenzyme active sites, including a model of the catalytic site of the zinc-bearing anthrax toxin lethal factor, a popular target for drug development. In particular, we find that the electrostatically embedded three-body (EE-3B) method is able to reproduce conventionally calculated bond-breaking energies in a series of pentacoordinate and hexacoordinate zinc-containing systems with an average absolute error (averaged over 25 cases) of only 0.98 kcal/mol.

Assessment and Validation of the Electrostatically Embedded Many-Body Expansion for Metal-Ligand Bonding.

The electrostatically embedded many-body method has been very successful for calculating cohesive energies and relative conformational energies of clusters, and here we extend it to calculate bond breaking energies for metal-ligand bonds in inorganic coordination chemistry. We find that, on average, the electrostatically embedded pairwise additive method is able to predict bond energies yielded by conventional full-system calculations done at the same level of theory to within 2.5 kcal/mol and that the electrostatically embedded three-body method consistently yields energies within 1.0 kcal/mol of the full-system calculations.

Evaluation of Density Functionals, SCC-DFTB, Neglect of Diatomic Differential Overlap (NDDO) Models and Molecular Mechanics Methods for Prolyl-Leucyl-Glycinamide (PLG) and Structural Derivatives.

Prolyl-leucyl-glycinamide (PLG) is a unique endogenous peptide that modulates dopamine receptor subtypes of the D(2) receptor family within the CNS. We seek to elucidate the structural basis and molecular mechanism by which PLG and its analogues modulate dopamine receptors, toward the development of new therapeutics to treat Parkinson's disease, tardive dyskinesia and schizophrenia. As a first step toward establishing a validated protocol for accurate computational modeling of PLG and associated peptidomimetic analogues, we evaluated the accuracy of density functional theory (DFT), wavefunction theory (WFT), and molecular mechanics (MM) calculations for PLG and for a library of structurally related small molecules. We first tested twelve local and nonlocal density functionals, Hartree-Fock (HF) theory, four "semiempirical" methods of the neglect of diatomic differential overlap (NDDO) type, and one self-consistent-charge nonorthogonal tight-binding (SCC-DFTB) method as implemented in two software suites, against coupled-cluster benchmark geometries for 4-methylthiazolidine, a small molecule that comprises key structural features present in our PLG analogue library. DFT and HF calculations were done with the MG3S augmented polarized triple-zeta basis set. We find that for 4-methylthiazolidine bond distances, DFT significantly outperforms NDDO, and both SCC-DFTB versions we evaluated perform worse than HF theory and are less accurate than 83% of the density functionals tested. The top five functionals for 4-methylthiazolidine were M05-2X, mPW1PW, B97-2, M06-2X, and PBEh, with mean unsigned errors (MUEs) in bond length of 0.0017, 0.0020, 0.0023, 0.0025 and 0.0027 Å, respectively. The widely used B3LYP functional ranked 11(th) out of twelve functionals evaluated, slightly below SCC-DFTB, and is significantly less accurate for 4-methylthiazolidine bond distances (MUE = 0.0095 Å) than the best local functional (M06-L, MUE = 0.0030 Å), which is far less computationally costly. Based on that initial analysis, we obtained new M05-2X benchmark geometric parameters for PLG and a library of eleven peptidomimetic derivatives, which we in turn used to examine the accuracy of thirty-four popular molecular mechanics (MM) force fields, four NDDO approaches, and SCC-DFTB for the full compound structures. Here, we found that ∼70% of the MM force fields tested superior to the best semiempirical and SCC-DFTB codings. Moreover, AMBER-type force fields proved most accurate among MM methods for this class of small-molecule peptidomimetics; the AMBER-type methods comprised eight out of the top ten molecular mechanics options we tested.

Design, synthesis and evaluation of analogs of initiation factor 4E (eIF4E) cap-binding antagonist Bn7-GMP.

Aberrant regulation of cap-dependent translation has been frequently observed in the development of cancer. Association of the cap-binding protein eIF4E with N(7)-methylated guanosine capped mRNA is the rate limiting step governing translation initiation; and therefore represents an attractive process for cancer drug discovery. Previously, replacement of the 7-Me group of the Me(7)-guanosine monophosphate with a benzyl group has been found to increase binding affinity to eIF4E. Recent X-ray crystallographic studies have revealed that the cap-binding pocket undergoes a unique structural change in order to accommodate the benzyl group. To explore the structure-activity relationships governing the affinity of N(7)-benzylated guanosine monophosphate (Bn(7)-GMP) for eIF4E, we virtually screened a library of 80 Bn(7)-GMP analogs utilizing CombiGlide as implemented in Schrodinger. A subset library of substituted Bn(7)-GMP analogs was synthesized and their dissociation constants (K(d)) were determined. Due to the poor correlation between docking/scoring results and experimental binding affinities, three-dimensional quantitative structure-activity relationship (3D-QSAR) calculations were performed. Two highly predictive and self-consistent CoMFA (comparative molecular field analysis) and CoMSIA (comparative molecular similarity indices analysis) models were derived and optimized. These models may be useful for the future design of eIF4E cap-binding antagonists.

Total synthesis and evaluation of C26-hydroxyepothilone D derivatives for photoaffinity labeling of beta-tubulin.

Three photoaffinity labeled derivatives of epothilone D were prepared by total synthesis, using efficient novel asymmetric synthesis methods for the preparation of two important synthetic building blocks. The key step for the asymmetric synthesis of (S,E)-3-(tert-butyldimethylsilyloxy)-4-methyl-5-(2-methylthiazol-4-yl)pent-4-enal involved a ketone reduction with (R)-Me-CBS-oxazaborolidine. For the synthesis of (5S)-5,7-di[(tert-butyldimethylsilyl)oxy]-4,4-dimethylheptan-3-one an asymmetric Noyori reduction of a beta-ketoester was employed. The C26 hydroxyepothilone D derivative was constructed following a well-established total synthesis strategy and the photoaffinity labels were attached to the C26 hydroxyl group. The photoaffinity analogues were tested in a tubulin assembly assay and for cytotoxicity against MCF-7 and HCT-116 cancer cell lines. The 3- and 4-azidobenzoic acid analogues were found to be as active as epothilone B in a tubulin assembly assay, but demonstrated significantly reduced cellular cytotoxicity compared to epothilone B. The benzophenone analogue was inactive in both assays. Docking and scoring studies were conducted that suggested that the azide analogues can bind to the epothilone binding site, but that the benzophenone analogue undergoes a sterically driven ligand rearrangement that interrupts all hydrogen bonding and therefore protein binding. Photoaffinity labeling studies with the 3-azidobenzoic acid derivative did not identify any covalently labeled peptide fragments, suggesting that the phenylazido side chain was predominantly solvent-exposed in the bound conformation.

Identification of novel non-hydroxamate anthrax toxin lethal factor inhibitors by topomeric searching, docking and scoring, and in vitro screening.

Anthrax is an infectious disease caused by Bacillus anthracis, a Gram-positive, rod-shaped, anaerobic bacterium. The lethal factor (LF) enzyme is secreted by B. anthracis as part of a tripartite exotoxin and is chiefly responsible for anthrax-related cytotoxicity. As LF can remain in the system long after antibiotics have eradicated B. anthracis from the body, the preferred therapeutic modality would be the administration of antibiotics together with an effective LF inhibitor. Although LF has garnered a great deal of attention as an attractive target for rational drug design, relatively few published inhibitors have demonstrated activity in cell-based assays and, to date, no LF inhibitor is available as a therapeutic or preventive agent. Here we present a novel in silico high-throughput virtual screening protocol that successfully identified 5 non-hydroxamic acid small molecules as new, preliminary LF inhibitor scaffolds with low micromolar inhibition against that target, resulting in a 12.8% experimental hit rate. This protocol screened approximately 35 million nonredundant compounds for potential activity against LF and comprised topomeric searching, docking and scoring, and drug-like filtering. Among these 5 hit compounds, none of which has previously been identified as a LF inhibitor, three exhibited experimental IC(50) values less than 100 microM. These three preliminary hits may potentially serve as scaffolds for lead optimization as well as templates for probe compounds to be used in mechanistic studies. Notably, our docking simulations predicted that these novel hits are likely to engage in critical ligand-receptor interactions with nearby residues in at least two of the three (S1', S1-S2, and S2') subsites in the LF substrate binding area. Further experimental characterization of these compounds is in process. We found that micromolar-level LF inhibition can be attained by compounds with non-hydroxamate zinc-binding groups that exhibit monodentate zinc chelation as long as key hydrophobic interactions with at least two LF subsites are retained.

Total synthesis and evaluation of 22-(3-azidobenzoyloxy)methyl epothilone C for photoaffinity labeling of beta-tubulin.

The total synthesis of 22-(3-azidobenzoyloxy)methyl epothilone C is described as a potential photoaffinity probe to elucidate the beta-tubulin binding site. A sequential Suzuki-aldol-Yamaguchi macrolactonization strategy was utilized employing a novel derivatized C1-C6 fragment. The C22-functionalized analog exhibited good activity in microtubule assembly assays, but cytotoxicity was significantly reduced. Molecular modeling simulations indicated that excessive steric bulk in the C22 position is accommodated by the large hydrophobic pocket of the binding site. Photoaffinity labeling studies were inconclusive suggesting non-specific labeling.

Energies, Geometries, and Charge Distributions of Zn Molecules, Clusters, and Biocenters from Coupled Cluster, Density Functional, and Neglect of Diatomic Differential Overlap Models.

We present benchmark databases of Zn-ligand bond distances, bond angles, dipole moments, and bond dissociation energies for Zn-containing small molecules and Zn coordination compounds with H, CH3, C2H5, NH3, O, OH, H2O, F, Cl, S, and SCH3 ligands. The test set also includes clusters with Zn-Zn bonds. In addition, we calculated dipole moments and binding energies for Zn centers in coordination environments taken from zinc metalloenzyme X-ray structures, representing both structural and catalytic zinc centers. The benchmark values are based on relativistic-core coupled cluster calculations. These benchmark calculations are used to test the predictions of four density functionals, namely B3LYP and the more recently developed M05-2X, M06, and M06-2X levels of theory, and six semiempirical methods, including neglect of diatomic differential overlap (NDDO) calculations incorporating the new PM3 parameter set for Zn called ZnB, developed by Brothers and co-workers, and the recent PM6 parametrization of Stewart. We found that the best DFT method to reproduce dipole moments and dissociation energies of our Zn compound database is M05-2X, which is consistent with a previous study employing a much smaller and less diverse database and a much larger set of density functionals. Here we show that M05-2X geometries and single-point coupled cluster calculations with M05-2X geometries can also be used as benchmarks for larger compounds, where coupled cluster optimization is impractical, and in particular we use this strategy to extend the geometry, binding energy, and dipole moment databases to additional molecules, and we extend the tests involving crystal-site coordination compounds to two additional proteins. We find that the most predictive NDDO methods for our training set are PM3 and MNDO/d. Notably, we also find large errors in B3LYP for the coordination compounds based on experimental X-ray geometries.

Total synthesis and evaluation of C25-benzyloxyepothilone C for tubulin assembly and cytotoxicity against MCF-7 breast cancer cells.

The total synthesis of C25-benzyloxy epothilone C is described. A sequential Suzuki-Aldol-Yamaguchi macrolactonization strategy was utilized employing a novel derivatized C8-C12 fragment. The C25-benzyloxy analog exhibited significantly reduced biological activity in microtubule assembly and cytotoxicity assays. Molecular modeling simulations indicated that excessive steric bulk in the C25 position may reduce activity by disrupting key hydrogen bonds that are crucial for epothilone binding to beta-tubulin.

Molecular Modeling of Geometries, Charge Distributions, and Binding Energies of Small, Drug-Like Molecules Containing Nitrogen Heterocycles and Exocyclic Amino Groups in the Gas Phase and Aqueous Solution.

We have tested a variety of approximate methods for modeling 30 systems containing mixtures of nitrogen heterocycles and exocyclic amines, each of which is studied with up to 31 methods in one or two phases (gaseous and aqueous). Fifteen of the systems are protonated, and 15 are not. We consider a data set consisting of geometric parameters, partial atomic charges, and water binding energies for the methotrexate fragments 2-(aminomethyl)pyrazine and 2,4-diaminopyrimidine, as well as their cationic forms 1H-2-(aminomethyl)pyrazine and 1H-2,4-diaminopyrimidine. We first evaluated the suitability of several density functionals with the 6-31+G(d,p) basis set to serve as a benchmark by comparing calculated molecular geometries to results obtained from coupled-cluster [CCSD/6-31+G(d,p)] wave function theory (WFT). We found that the M05-2X density functional can be used to obtain reliable geometries for our data set. To accurately model partial charges in our molecules, we elected to utilize the well-validated Charge Model 4 (CM4). In the process of establishing benchmark values, we consider gas-phase coupled cluster and density functional theory (DFT) calculations followed by aqueous-phase DFT calculations, where the effect of solvent is treated by the SM6 quantum mechanical implicit solvation model. The resulting benchmarks were used to test several widely available and economical semiempirical molecular orbital (SE-MO) methods and molecular mechanical (MM) force fields for their ability to accurately predict the partial charges, binding energies to a water molecule, and molecular geometries of representative fragments of methotrexate in the gaseous and aqueous phases, where effects of water were simulated by the SM5.4 and SM5.42 quantum mechanical implicit solvation models for SE-MO and explicit solvation used for MM. In addition, we substituted CM4 charges into the MM force fields tested to observe the effect of improved charge assignment on geometric and energetic modeling. The most accurate MM force fields (with or without CM4 charges substituted) were validated against gas-phase and aqueous-phase geometries and charge distributions of a larger set of 16 drug-like ligands, both neutral and cationic. This process showed that the Merck Molecular Force Field (MMFF94) with or without CM4 charges substituted, is, on average, the most accurate force field for geometries of molecules containing nitrogen heterocycles and exocyclic amino groups, both protonated and unprotonated. This force field was then applied to the complete methotrexate molecule, in an effort to systematically explore its accuracy for trends in geometries and charge distributions. The most accurate force fields for the binding energies of nitrogen heterocycles to a water molecule are OPLS2005 and AMBER.

Zn Coordination Chemistry: Development of Benchmark Suites for Geometries, Dipole Moments, and Bond Dissociation Energies and Their Use To Test and Validate Density Functionals and Molecular Orbital Theory.

We present nonrelativistic and relativistic benchmark databases (obtained by coupled cluster calculations) of 10 Zn-ligand bond distances, 8 dipole moments, and 12 bond dissociation energies in Zn coordination compounds with O, S, NH3, H2O, OH, SCH3, and H ligands. These are used to test the predictions of 39 density functionals, Hartree-Fock theory, and seven more approximate molecular orbital theories. In the nonrelativisitic case, the M05-2X, B97-2, and mPW1PW functionals emerge as the most accurate ones for this test data, with unitless balanced mean unsigned errors (BMUEs) of 0.33, 0.38, and 0.43, respectively. The best local functionals (i.e., functionals with no Hartree-Fock exchange) are M06-L and τ-HCTH with BMUEs of 0.54 and 0.60, respectively. The popular B3LYP functional has a BMUE of 0.51, only slightly better than the value of 0.54 for the best local functional, which is less expensive. Hartree-Fock theory itself has a BMUE of 1.22. The M05-2X functional has a mean unsigned error of 0.008 Å for bond lengths, 0.19 D for dipole moments, and 4.30 kcal/mol for bond energies. The X3LYP functional has a smaller mean unsigned error (0.007 Å) for bond lengths but has mean unsigned errors of 0.43 D for dipole moments and 5.6 kcal/mol for bond energies. The M06-2X functional has a smaller mean unsigned error (3.3 kcal/mol) for bond energies but has mean unsigned errors of 0.017 Å for bond lengths and 0.37 D for dipole moments. The best of the semiempirical molecular orbital theories are PM3 and PM6, with BMUEs of 1.96 and 2.02, respectively. The ten most accurate functionals from the nonrelativistic benchmark analysis are then tested in relativistic calculations against new benchmarks obtained with coupled-cluster calculations and a relativistic effective core potential, resulting in M05-2X (BMUE = 0.895), PW6B95 (BMUE = 0.90), and B97-2 (BMUE = 0.93) as the top three functionals. We find significant relativistic effects (∼0.01 Å in bond lengths, ∼0.2 D in dipole moments, and ∼4 kcal/mol in Zn-ligand bond energies) that cannot be neglected for accurate modeling, but the same density functionals that do well in all-electron nonrelativistic calculations do well with relativistic effective core potentials. Although most tests are carried out with augmented polarized triple-ζ basis sets, we also carried out some tests with an augmented polarized double-ζ basis set, and we found, on average, that with the smaller basis set DFT has no loss in accuracy for dipole moments and only ∼10% less accurate bond lengths.

Assessment of Density Functionals, Semiempirical Methods, and SCC-DFTB for Protonated Creatinine Geometries.

Creatinine concentrations in blood and urine can be used to detect renal insufficiencies and muscle diseases, but current chemical sensors cannot measure creatinine with sufficient selectivity and robustness because they lack a receptor that binds protonated creatinine (creatininium) selectively enough. As a first step toward identifying potential receptors for creatininium, we examine the accuracy of density functional theory (DFT) and wave function theory (WFT) calculations for creatininium cation geometries, evaluated against reference parameters from experiment. We tested twenty-one local and nonlocal density functionals, Hartree-Fock theory, four semiempirical molecular orbital (SEMO) methods of the neglect of differential overlap (NDO) type, and one self-consistent-field nonorthogonal tight-binding method (SCC-DFTB) as implemented in two closely related software packages. DFT and HF calculations were carried out with the MG3S augmented polarized triple-zeta basis set. We find that DFT significantly outperforms SEMO methods for our dataset, and both SCC-DFTB releases we tested (which gave almost identical results) were less accurate than 81% of the density functionals evaluated. The top five functionals for the creatininium structures we examined were MPW1B95, PBEh, mPW1PW, SVWN5, and B97-2, with MMUEs in bond length of 0.0126 Å, 0.0129 Å, 0.0133 Å, 0.0142 Å, and 0.0144 Å respectively, which indicates that all five functionals are similarly suitable for creatininium. The popular B3LYP functional has a MMUE of 0.0178 Å, which ranks it 16(th) overall. B3LYP also performs less favorably than the best local functional (SVWN5), which is less expensive.

A preliminary in silico lead series of 2-phthalimidinoglutaric acid analogues designed as MMP-3 inhibitors.

Matrix metalloproteinases (MMPs) have been the subject of intense research because of their roles in tumor metastasis and in the rise and spread of degenerative diseases such as osteo- and rheumatoid arthritis. A preliminary class of 140 druglike, small-molecule matrix metalloproteinase-3 inhibitors, intended as starting scaffolds for optimization and synthesis, has been designed in silico using a series of highly predictive three-dimensional quantitative structure-activity relationship models, including comparative molecular field analysis and comparative molecular similarity indices analysis, with docking and scoring. Thalidomide was chosen as the skeleton on which to base the new lead series, as it moderately inhibits MMP-3, is antiangiogenic, and lends itself easily to structural modifications. Most of the new compounds demonstrate medium to high predicted biological activity and good bioavailability as estimated by the octanol-water partition coefficient ClogP. Compound 102 in particular exhibits extremely favorable predicted activity against MMP-3; is moderately bioavailable; satisfies Lipinski's Rule of Five; and shows promise for further optimization, synthesis, and experimental evaluation as a potential adjunct anticancer or antirheumatic therapeutic.

Highly predictive CoMFA and CoMSIA models for two series of stromelysin-1 (MMP-3) inhibitors elucidate S1' and S1-S2' binding modes.

Three-dimensional quantitative structure-activity relationship models have been derived using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) for two training sets of arylsulfonyl isoquinoline-based and thazine/thiazepine-based matrix metalloproteinase inhibitors (MMPIs). The crystal structure of stromelysin-1 (MMP-3) was used to pinpoint areas on the ligands and receptors where steric and electrostatic effects (for CoMFA) and steric, electrostatic, hydrogen-bond donor, hydrogen-bond acceptor, and hydrophobic effects (for CoMSIA) correlate with an increase or decrease in experimental biological activity. The most predictive CoMFA and CoMSIA models were obtained using training-series subsets that sampled a wide range of activities, together with docking and scoring, inertial alignment, investigation of various partial charge formalisms, and manual adjustment of each compound within the active site. The models developed in this study are in agreement with experimentally observed MMP-3 structure-activity relationship data and offer new insights into binding modes involving the partly solvent-exposed S1-S2' subpocket and certain zinc-chelating groups.

Coordinative properties of highly fluorinated solvents with amino and ether groups.

Despite the widespread use of perfluorinated solvents with amino and ether groups in a variety of application fields, the coordinative properties of these compounds are poorly known. It is generally assumed that the electron withdrawing perfluorinated moieties render these functional groups rather inert, but little is known quantitatively about the extent of their inertness. This paper reports on the interactions between inorganic monocations and perfluorotripentylamine and 2H-perfluoro-5,8,11-trimethyl-3,6,9,12-tetraoxapentadecane, as determined with fluorous liquid-membrane cation-selective electrodes doped with tetrakis[3,5-bis(perfluorohexyl)phenyl]borate salts. The amine does not undergo measurable association with any ion tested, and its formal pK(a) is shown to be smaller than -0.5. This is consistent with the nearly planar structure of the amine at its nitrogen center, as obtained with density functional theory calculations. The tetraether interacts very weakly with Na(+) and Li(+). Assuming 1:1 stoichiometry, formal association constants were determined to be 2.3 and 1.5 M(-1), respectively. This disproves an earlier proposition that the Lewis base character in such compounds may be nonexistent. Due to the extremely low polarity of fluorous solvents and the resulting high extent of ion pair formation, a fluorophilic electrolyte salt with perfluoroalkyl substituents on both the cation and the anion had to be developed for these experiments. In its pure form, this first fluorophilic electrolyte salt is an ionic liquid with a glass transition temperature, T(g), of -18.5 degrees C. Interestingly, the molar conductivity of solutions of this salt increases very steeply in the high concentration range, making it a particularly effective electrolyte salt.